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Abstract 

Measurements and modeling of electron spin transport and dynamics are used to characterize 
hyperfine interactions in Fe/GaAs devices with n-GaAs channels. Ga and As nuclei are polarized 
by electrically injected electron spins, and the nuclear polarization is detected indirectly through 
the depolarization of electron spins in the hyperfine field. The dependence of the electron spin 
signal on injector bias and applied field direction is modeled by a coupled drift-diffusion equation, 
including effective fields from both the electronic and nuclear polarizations. This approach is used 
to determine the electron spin polarization independently of the assumptions made in standard 
transport measurements. The extreme sensitivity of the electron spin dynamics to the nuclear spin 
polarization also facilitates the electrical detection of nuclear magnetic resonance. 
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Recent observations of electrical spin injection and detection in ferromagnet - semi- 
conductor (FS) devices have been confirmed by their demonstrated sensitivity to electron 
spin precession.— There remain, however, several important issues which must be re- 
solved in order to interpret spin transport measurements. For example, in ra-GaAs doped 
near the metal-insulator transition, spin dynamics are influenced profoundly by hyperfine 
interactions.- Coupling between electron and nuclear spins leads to significant deviations 
from the expectations of simple drift-diffusion models.-^ More importantly, the quantita- 
tive determination of the electron spin polarization in a FS device using standard transport 
measurements is based on assumptions about boundary conditions and densities of states 
that are not directly verifiable. A detailed understanding of hyperfine interactions could 
allow for an independent measurement of the electron spin polarization. 

In this article we present measurements employing electrical generation and detection of 
dynamic nuclear polarization (DNP) in Fe/n-GaAs heterostructures. Non-equilibrium spin 
polarization is electrically injected into a GaAs channel, inducing DNP.- 1 ^^ We show that 
the resulting dynamics of the combined electron-nuclear spin system can be described in a 
completely self-consistent manner. This approach allows for a measurement of the electron 
spin polarization that is independent of assumptions about either boundary conditions at 
the Fe/GaAs interface or the electronic densities of states in the ferromagnet or the semi- 
conductor. Finally, we show that the spin-polarized electrons in GaAs in combination with 
the ferromagnetic detector can be used as a sensitive probe of nuclear magnetic resonance. 

A schematic diagram of one of our devices is shown in Fig. [T^a).-^ 1 ^ The heterostructures 
consist of a 2.5 /im thick Si-doped n-GaAs (n = 5 x 10 16 cm~ 3 ) channel and 5 nm thick 
Fe electrodes that are deposited epitaxially on the GaAs (100) surface. A Schottky tunnel 
barrier is formed by a n —>■ n + (~5x 10 18 cm~ 3 ) GaAs transition layer— between the Fe and 
the n-type channel. The Fe injection and detection contacts are 7 x 50 /im 2 and 4 x 50 fim 2 
respectively, with a center to center gap of 9.5 /im. DNP is observed in all heterostructures 
of this general design below 80 K. At higher temperatures, spin-polarized electrons are no 
longer bound to the donor sites, thus suppressing the hyperfine coupling.- The discussion in 
this paper is based on data obtained at T = 60 K. 

Measurements are performed in the non-local 2 - 1 ^ 3 - or 3-terminalA 1 A configurations shown in 
Fig. []Ja). Under forward bias, electrons tunnel from GaAs into the Fe injector (contact b). 
Majority spins aligned in the x direction accumulate at the injector interface^ and diffuse 



2 








■(b) 



<L -40 - 



t i 



> 



-80 - 



T = 60 K | 

V- ■ = 0.35 V I. 
inj i 



Vu 



-300 







300 



B x (Oe) 



FIG. 1: (Color online) (a) Schematic of the lateral spin transport device. I establishes a forward 
bias at contact b. The spin polarization in the GaAs channel at contact b or c is determined 
by measuring the 3-terminal voltage V%t or non-local voltage Vnl relative to reference contact d 
(drawing not to scale), (b) Vnl for longitudinal field B x swept from large positive field (black line) 
and from large negative field (red dotted line). 

to the non-local detector (contact c). Contact d, located 290 fim away, serves as a voltage 
reference. A longitudinal magnetic field B x is swept along the Fe easy axis ([110] direction, 
labeled x in Fig. Ufa)) at a rate of 0.025 Oe/s in order to ensure that the nuclear spin 
polarization is in equilibrium. A small static transverse field of B z = 18 Oe is applied in the 
z direction. This field does not perturb the magnetization of the contacts. 

We first consider non-local spin-valve measurements. We observe jumps in the non-local 
voltage Vnl when the injector and detector magnetizations switch between parallel and anti- 
parallel states, as shown in Fig. [T](b) after subtraction of a spin-independent background. 
The magnitude of the voltage jump — V\i is proportional to PcaAs, where PcaAs = 
(n-f —ni)/ (n^+ni) and n-f and are the densities of up and down spins respectively. Under 
forward bias PcaAs increases with increasing Vi n j, as shown in Fig. [2J 

The other important aspects the data of Fig. [TJb) are the depolarization dip (Vnl ~ 
Vn)/2 and re-polarization peak (Vnl — Vf\) & t l° w field (\B X \ < 150 Oe). These features 
are shown on an expanded scale for several interfacial bias voltages in Figs. [3]^a) and (b). 
A three-terminal measurement^ is sensitive only to accumulation of spins parallel to the 
magnetization of the injector itself and hence is a direct probe of spin precession in the 
semiconductor. The fact that the field dependence of the two types of measurements are 
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very similar indicates that the low-field signatures are due entirely to precession. It can 
be easily verified, however, that the effective field leading to the precessional dynamics is 
much larger than the applied field. This reflects the existence of a hyperfine field due to 
dynamically polarized nuclei, which has the following form:— 

B rt =„J B t b ' §) : S(S+KS) (i) 

(B + KSr + (Bf, 

where b n and b e , which are both negative in GaAs,— represent effective fields due to the 
polarized nuclei and electrons, S is the average electron spin (\S\ = 1/2 for PcaAs = 100%), 
Bi is the local dipolar field experienced by the nuclei, and £ parameterizes the assisting 
processes which allow energy to be conserved in mutual spin flips between electrons and 
nuclei.— \Bn\ can be as large as several Tesla in our samples. 

From Eq. [U we determine that at large B x , Bn is essentially anti- aligned with S , as 
shown in Fig. [3(c), and S remains polarized, with negligible dephasing from precession. 
However, when B x becomes comparable to B Z1 B^ rotates towards the shown in 

Fig. [3(d). S precesses around B^, z resulting in the observed depolarization. At a smaller 
field B r , indicated by the arrow in Fig. [3(a), the electron spin system becomes re-polarized, 
indicating that the nuclear field has been suppressed. This phenomenon has previously been 
associated with the cancellation of the magnetic field along the spin injection axes B x by 
the electronic exchange field b e S acting on the nuclei.— When B x is swept from the opposite 
direction the spin polarization S and hence the Knight field are opposite in sign, and the 
sign of B r reverses, as shown in Fig. [1(b). 

In optical orientation experiments, the argument of the previous paragraph has been 
used to determine b e by setting B r = b e S, where S is the optically generated spin 
accumulation.— 1 ^ In our experiment, we can enhance S by increasing the injector bias. As 
shown in Fig. [2]( inset), however, B r as measured by either the non-local or three-terminal 
methods clearly decreases with increasing bias, in direct contradiction of the bias depen- 
dence of Vj-f — V|j shown in Fig. [2j This anomaly, which precludes a simple identification 
of the Knight field, is due to the fact that the electron spin polarization at small fields is 
reduced by precession. This fact, which was ignored in the interpretation of the early optical 
orientation measurements, has a profound influence on the electron spin dynamics at very 
low fields. 

To address this shortcoming as well as other aspects of our experiments, including the 
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FIG. 2: (Color online) Vj-f — Vfj, (black open squares, left axis) and modeled PgoAsj the spin 
polarization averaged over the area of the non-local detector, i.e. under contact c, (red circles, 
right axis) as a function of injector bias. Inset: re-polarization field B r , as defined in the text, for 
non-local (blue solid squares) and 3-terminal (black open circles) measurements as a function of 
injector bias. 

spatial inhomogeneity in both the electron and nuclear spin polarizations, we consider a 
more complete model based on the drift-diffusion formalism that has been applied widely to 
spin transport problems.-^^^^ The dynamics of the injected spins in the GaAs channel 
are given by the following rate equation: 



dS 
~dt 



ds a 2 s s =, 

-v d — + D-— Q L x S, 



(2) 



' dy dy 2 r. 

where v,i, D, and t s are the drift velocity, diffusion constant, and spin lifetime; Ql — 
gfiBBtott/h is the Larmor frequency, g = —0.44 is the electron g-factor in GaAs, \xb is the 
Bohr magneton, and B to t = B + Bn is the sum of the external applied field and the hyperfine 
field defined in Eq. [TJ In order to model the experimental data, we solve Eq. [2] using the 
Crank-Nicholson method^ with a one-dimensional spatial grid. A constant spin generation 
rate So is introduced in each cell beneath the injector. S and B^ are calculated at each 
position for each time step, thus leading to a self-consistent steady-state solution (typically 
after 50 ns).— The non-local and 3-terminal curves as a function of B x are determined from 
the spatial average of S x at the detector and injector respectively. The parameters Vd and D 
are determined independently from Hall and resistivity measurements. In the simulations, 
we set Vd — outside of the charge current path, i.e. between the injector and detector. 
b n ,b e , y/^Bi, and r s are obtained from fits of a field sweep at one bias current and are then 
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FIG. 3: (Color online) 3-terminal voltage V%t (a) and non-local voltage Vnl (b) as a function of 
applied field B x in the measurement configuration shown in Fig. [2(a) f° r different injector biases 
Vmj. Open circles are experimental data for B x swept from large negative field at T = 60 K. 
Solid lines are modeling results as described in the text. Schematics showing the direction of the 
hyperfine nuclear field for large B x (c), and for B x comparable to the transverse field B z (d). 



kept fixed while fitting the data for other bias currents. Only S varies with bias, but it is 
kept the same for simulations of the same bias. As indicated by the solid curves in Figs. [3(a) 
and (b), the modeling clearly reproduces the measured curves and the re-polarization fields.— 



From the modeling we determine that b n = —53 kOe and b e = —50 Oe, which are, as 
expected, smaller than the theoretical ideal values of b n = —170 kOe and b e = —170 Oe 
(calculated by Paget et alM for a closed electron-nuclear spin system and in which the donor 
sites are always occupied by spin-polarized electrons) but larger by a factor of 1.3 to 4 then 
previously measured values for p-type GaAs^ and Ge-doped GaAs.— We find y/^Bi = 40 Oe. 
£ incorporates sample specific processes that assist nuclear spin polarization and is given by 
£ = Tpoi/T^B/Bi) 2 . 10 Using T pol » 10 s,-^ B { = 1.45 Oe,i£ and a measured T x = 40 s at 
B = 100 Oe, we find y/^Bi = 50 Oe, which is comparable to our measured value. 

The complete consistency among the different curves in Fig. [3] given a single set of param- 
eters allows us to use this approach to measure the bias dependence of the spin polarization 
PGaAs = 2(|^|) averaged over the non-local detector. The results extracted from the mod- 
eling, shown in Fig. El are in excellent qualitative agreement with the bias dependence of 
Vj-f — V\i- It is extremely important to note that this measurement of PGaAs, which is based 
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FIG. 4: (Color online) Non-local voltage field dependence: measured curves (red circles)and 
modeled curves with Knight field (black line) and without Knight field (b e = 0, green dashed line) 
for magnetic field swept along the z axes (Hanle) (a) and along polar angle 6 = 17° from the z 
axes (oblique Hanle) (b). Modeled curves are calculated with the same parameters used to fit data 
in Fig. El (b) Inset: closeup of the oblique Hanle curves, (c) The three-terminal voltage VJjt as 
a function of frequency of applied ac magnetic field. The observed resonances for each isotope 
are indicated in the legend. Transitions at the sums of the fundamental frequencies of different 
isotopes are also observed. 

on the Knight field, is independent of any of the usual assumptions underlying spin trans- 
port measurements. It is therefore of particular interest to compare, quantitatively, with 
the polarization as inferred from the non-local voltage using standard arguments. Following 
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PcaAs can be estimated from — V\i = 2r]PF e PGaAsE / /3e, where rj is the interface 
spin transmission efficiency, Pp e = is the spin polarization of Fe, and Ef is the fermi 
energy of GaAs, which is assumed to be the that of a Pauli metal with an effective mass 
m* = 0.07m e . From the measured — in Fig. [2] we have calculated PcaAs (for rj = 0.5, 
which is expected based on spin-LED measurements^ and a measured carrier density of 
4.3 x 10 16 cm~ 3 at T = 60 K) to be between 0.2% and 6% for the injector bias range 
shown. This is in good agreement with the values obtained from the analysis of the coupled 
electron-nuclear spin dynamics shown in Fig. [2j It is not obvious a priori that this should 
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be the case, since the spin-polarized electrons responsible for the Knight shift are bound on 
donors while those probed by the non-local measurement are at the Fermi level. Among 
the important implications of this result is that the magnitude of the density of states in 
the semiconductor near the Fe/GaAs interface is not grossly different from that given by 
the naive Pauli model. We have also implemented this analysis on a separate Fe/GaAs het- 
erostructure with the same channel doping but with Vj-j — an order of magnitude lower. 
Fits with the same parameters as those used to model data in Fig. [3] yielded PcaAs that are 
in agreement with the significantly lower spin signal. 

The extraordinary sensitivity of the electron spin dynamics to hyperfine effects in the FS 
system was initially observed in the distortion of the Hanle effect, in which electron spin 
polarization is suppressed by precession in a transverse magnetic field.- Fig. H] shows V^l for 
transverse field swept along the z-axis and at a polar angle of 17° from the z-axis (the static 
field in Fig. [T](a) is eliminated). The salient features of these data are reproduced by the 
model introduced above with exactly the same parameters used to fit the data of Fig. [3l Of 
particular interest are the narrowing of the Hanle curve in Fig. IU(a) and the weak splitting of 
the zero-field peak shown in the inset of Fig. H](b). Both are due to the effect of the Knight 
field. In its absence, the Hanle curves in Fig. Ufa) and (b) would show the form indicated 
by the dashed curves. We believe that the model's overestimation of the magnitude of the 
high-field satellite peak in Fig. IH(b) is due to the assumptions of purely one-dimensional 
diffusion as well as a uniform current density across the injection contact. 

A final indication of the strong coupling of the electron and nuclear spin systems is 
provided by the detection of nuclear magnetic resonance (NMR) when an ac magnetic field 
is applied by a small coil placed over the sample. This is illustrated in Fig. H^c), in which 
V%t is shown as a function of frequency in a constant oblique magnetic field at T = 30 K. 
When the nuclei are off-resonance, the electron spin polarization is suppressed by precession 
around B^. At a resonance, the nuclei are partially thermalized, is reduced, and the 
electron spin polarization is restored. At low fields, this approach can be used to observe 
all of the possible transitions in the 1 = 3/2 manifold for each of the three isotopes ( 69 Ga, 
71 Ga, and 75 As) present in the sample. The observation of higher order (Am = ±2, ±3) 
transitions is likely due to dipolar interactions (which also lead to the "sum transitions" 
from two different isotopes) or quadrupolar coupling. 23 

In conclusion, we have demonstrated that the coupled electron-nuclear spin dynamics in 



S 



Fe/n-GaAs spin transport devices can be understood quantitatively, providing an indepen- 
dent means of determining the electron spin polarization. This work was supported by NSF 
under DMR 0804244, the ONR MURI program, the NSF NNIN program, and the Japan 
Society for the Promotion of Science (T.K.). 



* On leave from Tokyo Institute of Technology, Yokohama, Kanagawa 226-8503, Japan. Current 
address: Corporate R&D center, Toshiba Corp., Kawasaki, Japan 

1 X. Lou, C. Adelmann, M. Furis, S.A. Crooker, C.J. Palmstr0m, and P.A. Crowell, Phys. Rev. 
Lett. 96, 176603 (2006). 

2 X. Lou, C. Adelmann, S.A. Crooker, E.S. Garlid, J. Zhang, K.S.M. Reddy, S.D. Flexner, C.J. 
Palmstr0m, and P.A. Crowell, Nat. Phys. 3, 197 (2007). 

3 I. Appelbaum, B. Huang, and D. Monsma, Nature (London) 447, 295 (2007). 

M. Tran, H. Jaffres, C. Deranlot, J. M. George, A. Fert, A. Miard, and A. Lemaitre, Phys. Rev 
Lett. 102, 036601 (2009). 

5 C. Awo-Affouda, O.M.J, van 't Erve, G. Kioseoglou, AT. Hanbicki, M. Holub, C.H. Li, and 
B.T. Jonker, Appl. Phys. Lett. 94, 102511 (2009) 

6 M. Ciorga, A. Einwanger, U. Wurstbauer, D. Schuh, W. Wegscheider, and D. Weiss, Phys. Rev. 
B 79, 165321 (2009). 

7 G. Salis, A. Fuhrer, and S. F. Alvarado, arxiv:cond-mat/0908.0273 

8 Optical Orientation, edited by F. Meier and B.P. Zakharchenya (North-Holland, New York, 
1984). 

9 A.W. Overhauser, Phys. Rev. 92, 411 (1953). 

10 D. Paget, G. Lampel, B. Sapoval, and V.I. Safarov, Phys. Rev. B 15, 5780 (1977). 

11 S.A. Crooker et. al, Science 309, 2191-2195 (2005). 

12 AT. Hanbicki et. al, Appl. Phys. Lett. 82, 4092-4094 (2003). 

13 M. Johnson and R.H. Silsbee, Phys. Rev. Lett. 55,1790 (1985hPhys. Rev. B 37, 5326 (1988). 

14 We adopt a notation similar to that used by Paget et al. in Ref. 10|, except that we have excluded 
additional numerical prefactors for b n and b e (the leakage factor / and effective donor occupancy 
r) which cannot be determined independently. 

15 B.P. Zakharchenya, V.K. Kalevich, V.D. Kul'kov, and V.G. Fleisher, Fiz. Tverd. Tela 



9 



(Leningrad) 23, 1387-1394 (1981). 

16 F. J. Jedema, H. B. Heersche, A. T. Filip, J. J. A. Baselmans, and B. J. van Wees, Nature 
(London) 416, 713 (2002) 

17 W.H. Press, S. A. Teukolsky, W.T. Vetterling, and B.P. Flannery Numerical Reciepes in C, The 
Art of Scientific Computing 2 nd ed., (Cambridge University Press, 1992). 

18 This approach ignores the dynamics of the nuclear spin system, for which Eq. Q] is true only in 
equilibrium. It will nevertheless arrive at a correct self-consistent steady-state solution for the 
coupled electron-nuclear spin dynamics, provided that the solution at equilibrium is unique. 

19 In fitting the data, the error in So is estimated to be about 20%. The error in nuclear field 
parameters is less than 20%. 

20 M.I. Dyakonov and V.I. Perel, in Optical Orientation, edited by F, Meier and B.P. Zakharchenya 
(North-Holland, New York, 1984) 

21 R.J. Soulen Jr. et. al, Science 282, 85-88 (1998). 

22 C. Adelmann, X. Lou, J. Strand, C.J. Palmstrom, and P.A. Crowell, Phys. Rev. B 71,121301(R) 
(2005). 

23 J. Strand, X. Lou, C. Adelmann, B.D. Schultz, A.F. Isakovic, C.J. Palmstr0m, and P.A. Crowell, 
Phys. Rev. B 72 155308 (2005). 



10 



